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We reveal the hydrogen isotope effect of three chemical reactions, i.e, the reflection, the absorption and the 
penetration ratios, by classical molecular dynamics simulation with a modified Brenner's reactive empirical 
bond order (REBO) potential potential. We find that the reflection by n— electron does not depend on the mass 
of the incident isotope, but the peak of the reflection by nuclear moves to higher side of incident energy. In 
addition to the reflection, we also find that the absorption ratio in the positive z side of the graphene becomes 
larger, as the mass of the incident isotope becomes larger. On the other hand, the absorption ratio in the negative 
z side of the graphene becomes smaller. Last, it is found that the penetration ratio does not depend on the mass 
of the incident isotope because the graphene potential is not affected by the mass. 
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1 Introduction 

Plasma-carbon interaction yields small hydrocarbon molecules on divertor region of a nuclear fusion deviceJU 
El El H |21 . Diffusing from divertor region to core plasma region of fusion device, generated hydrocarbon takes 
energy from the core plasma. Reduction of hydrocarbon diffusing from divertor is the main aim of studies in 
plasma-carbon research. To achieve the aim, researches with computer simulation have been being done |6] [7] [8] 
[9]. However, the creation mechanism of the hydrocarbons has not been elucidated yet. 

We, therefore, as the first step to clarify the creation mechanism, investigated, by computer simulation, colli- 
sion process of hydrogen atoms and one graphene sheet, which is regarded as one of basic processes of complex 
plasma-carbon interaction in the previous works[7. 8|. From the previous works in which an incident hydrogen 
kinetic energy E\ is less than 100 eV to compare with experiments, it was found that an hydrogen-absorption 
ratio of one graphene sheet depends on the incident hydrogen energy, and that the collision mechanism between 
a graphene and a hydrogen can be classified into three types of processes: absorption process, reflection process, 
and penetration process (see Fig. [TJ a ))- Moreover, it was also found that when hydrogen atom is absorbed by 
graphene, the nearest carbon atom overhangs from the graphene which we called "overhang structure". 

Based on the above results, as the second step, simulation model were extendedQ from a single graphene 
sheet to multilayer graphene sheets, which is a more realistic sputtering process of graphene sheets and hydrogen 
atoms than the previous workQ. From the second workjg), we found the following fact: breaking the covalent 
bonds between carbon atoms by hydrogen does not play an important role during destruction process of graphene 
structure, but momentum transfer from incident hydrogen to graphene causes to destroy graphene structure. 
Moreover, it was found[(9), that almost all fragments of graphene sheets form chain-shaped molecules, and that 
yielded hydrocarbon molecules are composed of carbon chain and single hydrogen-atom. 

In the present paper, we investigate hydrogen isotope effect for collision process of a single hydrogen isotope 
and one graphene sheet. Information of dependence of the chemical reaction on a type of isotope is necessary 
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to realize plasma confinement nuclear fusion system. In the present simulation, we change only the mass of 
the injected isotope, without changing the interaction potential. We used 'classical' molecular dynamics (CMD) 
algorithm with modified Brenner's reactive empirical bond order (REBO) potential which we proposed to deal 
with chemical reaction between hydrogen and graphene in the previous simulations 171 [8l l9l [TUl . 
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Fig. 1 Incident energy dependence of the absorption, the reflection and the penetration ratios. Three types of injected atoms 
are hydrogen (a), deuterium (b) and tritium (c). Open triangles, open circles, and open squares denote denote the absorption, 
the reflection and the penetration ratios, respectively. Dash-dotted lines, solid lines, and short-dashed lines are drawn as the 
guide for eyes for the absorption, the reflection and the penetration ratios, respectively. 



2 Simulation Method and Model 

We adapt CMD simulation with the NVE condition, in which the number of particles, volume and total energy 
are conserved. The second order symplectic integration ifTTl is used to solve the time evolution of the equation of 
motion. The time step is 5 x 10~ 18 s. The modified Brenner's reactive empirical bond order (REBO) potential|10j 
has the following form: 

u = E [^](^-)-^(W>^ B }>{0)^](^)" , 

i.j>i 

(1) 
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where ry is the distance between the i-th and the j-th atoms. The functions VSi and VjA, represent repulsion 
and attraction, respectively. The function 6^ generates multi-body force. (See details of the modified Brenner's 
REBO potential in Ref.JU.) 

In order to investigate the difference of the isotopes, i.e., hydrogen (H), deuterium (D) or tritium (T), we 
clarify the mass dependence of the injected isotope. The value of the mass for H, D, or T is shown in Table Q] 
The potential function is not changed for each isotope. 

Simulation model is shown in Fig. [2] We inject the hydrogen isotope into the graphene composed of 160 
carbon atoms. The center of mass of the graphene is set to the origin of coordinates. The surface of the graphene 
is parallel to the x-y plane. The size of the graphene is 2.13 nm x 1.97 nm. The graphene has no lattice defects 
and no crystal edges due to periodic boundary condition toward x and y directions. The structure of the graphene 
is used to the ideal lattice of graphene. Each velocity of carbon atoms of the graphene is zero in the initial state, 
that is, the initial temperature of the graphene is set to zero Kelvin. 

The hydrogen isotope is injected parallel to the z axis from z = 4 A. We repeat 200 simulations where the x 
and y coordinates of injection points are set at random. As a result, we obtain three chemical reaction ratios for 
H, D, or T by counting each a reaction. 



3 Results and Discussions 

We observed three kinds of reactions between the single hydrogen isotope atom and the graphene by CMD 
simulation, which are absorption reaction, reflection reaction and penetration reaction (see Fig. [T). We found 
the following differences of the reflection and the absorption ratios among three isotopes. On the other hand, the 
penetration ratio has almost the same E\ dependence. 



3 . 1 Reflection ratio 

From the previous workHfl. it was found that two kinds of repulsive force work between the incident atom and 
the graphene. One is derived by the 7r— electron over the graphene and the other is done by nuclear of carbon. 

As the result of the present simulation, the E\ dependence has the following properties. In the case of E\ < 
0.5 eV, the reflection ratio is almost one for all isotopes. This behavior is explained by the fact that the reflection 
in this energy region is derived by the repulsive force of tt— electrons over graphene surface|8j), which does not 
depend on the mass of the incident isotope. As E\ is getting larger than 0.5 eV, the reflection ratios are decreasing, 
and then it is increasing by the nuclear repulsive force of the carbon atom. Around E\ ~ 15 eV, they have the 
peak (see Figs. [T]and[3]l. Then, without the mass dependence, they decrease again and approach to zero (Fig. 
[3]), because the penetration reaction appears in the energy region that E\ > 15 eV. The details of the penetration 
reaction will appear in ^3.31 

By comparison with three isotopes, it is found that the peak energy of the reflection ratio becomes larger as 
the mass is getting larger, but the peak height becomes smaller (see Fig. 01. This behavior can be explained 
by the reflection mechanism in the previous work|8|, where the incident energy E\ has the following necessary 
condition for the reflection reaction by the nuclear of carbon: 

Ei > EretOm) = 0.84 eV , (2) 

?77h nic + m 

where the isotope mass m dependence is modified from the previous equation[8'| in order to cover three isotopes. 
The hydrogen mass toh is 1.00794 u, and the carbon mass mc is 12.00000 u. The value of E le f(m) is given in 
Table [T|and Fig. |3]for all isotopes. 

From Eq. (O or Table Q] it is found that, as the mass of the incident atom becomes larger, it needs higher 
incident energy for the isotope to be reflected by the graphene. The reflection reaction by the nuclear repulsion 
occurs at E\ ~ E m f, and then it increases monotonically until the penetration reaction becomes dominant among 
three reactions, as shown in Fig. |3]Therefore, in the case that the starting energy of the reflection reaction by the 
nuclear E Te f is smaller, the reflection ratio can reach a higher value at the peak energy. According to the above 
mechanism, the peak energy of the reflection ratio becomes larger and the peak height becomes smaller, as the 
mass of the isotopes is getting heavier, as shown in Fig. [3] 
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3.2 Absorption ratio 

The absorption ratio has two peaks at E\ ~ 5eV and 24 eV. One peak denotes the overhang state in the positive 
z side of the graphite, and the other peak is the overhang state in the negative z side[8 j. From Fig. [T| the height 
of the second peak, which is near 24 eV, becomes smaller, as the mass of the incident isotope is increasing. On 
the other hand, the first peak of the absorption ratio, which is around 5 eV, becomes large. The origin of the mass 
dependence of the absorption ratio is the same one as that of the reflection ratio in fJXT] The velocity becomes 
slower as the mass becomes large. Therefore, it becomes easier for the graphene to trap the isotope in the positive 
z side. 



3.3 Penetration ratio 

From Fig. Q] it is found that the penetration ratio does not depend on the mass of the incident isotope, i.e., H, D 
and T. The incident atom must overcome the graphene potential to penetrate the graphene sheet. The graphene 
potential does not depend on the mass of the incident atom. Therefore, the penetration ratio does not depend on 
the mass of the isotope. 



4 Summary 

We reveal the hydrogen isotope effect of three chemical reactions, i.e, the reflection, the absorption and the 
penetration reactions, by CMD simulation with the modified Brenner's REBO potential. From the previous 
workj8], the reflection process is divided into two processes, that is, reflection by it— electron and by nuclear. In 
the present work, we find that the reflection by it— electron does not depend on the mass of the incident isotope, 
but the peak of the reflection by nuclear moves to higher side of E\. In addition to the reflection, we also find 
that the absorption ratio in the positive z side of the graphene becomes larger, as the mass of the incident isotope 
becomes larger. On the other hand, the absorption ratio in the negative z side of the graphene becomes smaller. 
Last, it is found that the penetration ratio does not depend on the mass of the incident isotope because the graphene 
potential is not affected by the mass. 




Fig. 2 Simulation model. There are 160 carbon atoms 
and an injected hydrogen isotope. The length of the co- 
valent bond between carbons is 0. 142nm at temperature 
T — OK. The periodic boundary condition is used in x 
and y directions. A hydrogen atom is injected parallel to 
the z axis from z =0.4nm. 
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(^(Ed)(Et) Incident energy: £, [eV] 

Fig. 3 The isotope mass m dependence of the reflec- 
tion ratio by the nuclear repulsion. Filled circles, filled 
squares, and filled triangles denote the reflection ratios 
for hydrogen, deuterium and tritium, respectively. These 
plotted data are picked up from Fig. [TJ The estimated 
minimum nuclear-reflection-energies Eu, -Ed and Ej are 
given in Table [T] Dashed line, solid line and dash-dotted 
line are drawn as the guide for eyes. 
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Table 1 The mass and the reflection threshold energy by nulcear repulsive force E Te { for hydrogen, deuterium and tritium. 
The reflection threshold energy E Ie f depends on the mass of the isotope as Eq. Q2). We use the unified atomic mass u as the 
unit of mass. 





Hydrogen 


Deuterium 


Tritium 


m [u] 


1.00794 (= m H ) 


2.01410(=m D ) 


3.01605(=m T ) 


E Ici (m) [eV] 


0.84(= Eh) 


1.56(=£ D ) 


2.18(= E T ) 
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